Anisotropic flow in y/s = 2.76 TeV Pb+Pb collisions at the LHC 
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The results on elliptic flow in y/s — 2.76 TeV Pb+Pb collisions at the Large Hadron Collider 
(LHC) reported by the ALICE collaboration are remarkably similar to those for y/s = 200 GeV 
gold-gold collisions at the Relativistic Heavy Ion Collider (RHIC). This result is surprising, given 
the expected longer lifetime of the system at the higher collision energies. We show that it is 
nevertheless consistent with 3+1 dimensional viscous event-by-event hydrodynamic calculations, and 
demonstrate that elliptic flow at both RHIC and LHC is built up mostly within the first ~ 5fm/c of 
the evolution. We conclude that an "almost perfect liquid" is produced in heavy-ion collisions at the 
LHC. Furthermore, we present predictions for triangular flow as a function of transverse momentum 
for different centralities. 



The LHC era has barely begun, yet it is already pro- 
ducing significant physics results. In particular, the 
ALICE collaboration has demonstrated that the QGP 
system produced at y/s = 2.76 TeV at the LHC is in 
many ways similar to the system produced by RHIC at 
the much lower y/s = 200 GeV. Specifically, the elliptic 
flow, measured by the ALICE collaboration [l[ is surpris- 
ingly close to that measured by the STAR collaboration. 
Naively, one would have expected that the higher pres- 
sure and the longer lifetime of the QGP would make the 
effect of spatial anisotropy much greater at the LHC. 
To determine whether the LHC data is showing truly 
unexpected feature requires thorough analysis of hydro- 
dynamics with energy and entropy density appropriate 
for the LHC. In this work, we use a 3+1D viscous hy- 
drodynamic simulation model [H, Q to show that this is 
actually a natural consequence of the self- quenching of 
elliptic flow when the spatial eccentricity comes down 
below 0.1. 

As has been extensively reviewed 0, HJ, the elliptic 
flow, defined by the coefficient of cos 2<j) in the momentum 
distribution as follows 

dN__ 1 dN t 1 + 2 v cos2( /, + ) (1) 
dyd 2 pT 2-7T dyprdpr 

has been one of the most im 
quark-gluon plasma (QGP) 



jortant evidence of the 
A significant amount 
of theoretical work has been carried out by various groups 
for RHIC and some for LHC [IM3- Recently, we 
have emphasized the role of event-by-event fluctuations 
of the initial condition and the finite viscosity in a full 
3+1D hydrodynamic calculation Q in understanding the 
details of elliptic flow measured at RHIC. In this study 
we demonstrate the importance of both at LHC energies. 
In this work, we use a variant of the Israel-Stewart 



formalism jl4H18l ] derived in [H 
tensor is decomposed as 



] , where the stress-energy 
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where 



= (e + V)u»u v - Vg 



(2) 



(3) 



is the ideal fluid part with flow velocity w M . The local en- 
ergy density and local pressure are e and V ', respectively, 
and g 1 *" = diag(l, —1, —1, —1) is the metric tensor. The 
flow velocity is defined as the time-like eigenvector of T-^ 

T%u v = eu» (4) 

with the normalization u v u v = 1. The pressure is deter- 
mined by the equation of state 

V = V(e,p B ). (5) 

The evolution equations are d^T^ = and 



TV 



-w^(d a u a ), 

(6) 



where 



= 77 V 



(7) 



with A^"- 
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the local 3-metric, V = A» v d v 



G 



cd 



(8) 
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the local space derivative, and r\ the shear viscosity. 

In the r, rj s coordinate system used in this work, these 
equations can be re- written as hyperbolic equations with 
sources 

daTtf = -d a W ab + F b 

and 

d a (u a W cd ) = -(l/r 7r )(iy cd - S cd ) 

where F b and G cd contain terms introduced by the coor- 
dinate change from t, z to r, rj s as well as those introduced 
by the projections in Eq.©. 

To solve this system of equations, the 3+1D relativis- 
tic hydrodynamic simulation MUSIC 0, [3 is used. This 
approach utilizes the Kurganov-Tadmor (KT) scheme 
[20I . |2~J ] , together with Heun's method to solve resulting 
ordinary differential equations. 

The initialization of the energy density is done using 
the Glauber model (see [22| and references therein): Be- 
fore the collision the density distribution of the two nuclei 
is described by a Woods-Saxon parametrization 

Po 



Pa(t) 



1 + cxp[(r - R)/d] 



(10) 
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with R — 6.62 fm and d = 0.546 fm for Pb nuclei. The 
normalization factor p is set to fulfil J d 3 rpA(r) = A. 
The relevant quantity for the following considerations is 
the nuclear thickness function 



Ta(x,v) 



dz p A {x,y,z) 



(11) 



where r = \fa 



The opacity of the nucleus 



is obtained by multiplying the thickness function with 
the total inelastic cross-section a of a nucleus- nucleus 
collision. 

The initial energy density distribution in the trans- 
verse plane is scaled with the number of wounded nu- 
cleons JiwN- The initial energy density at the center 
is 175GeV/fm 3 in the ideal case, and 158GeV/fm for 
t]js = 0.08: the ratio of shear viscosity to entropy den- 
sity. For the event-by-event simulation, the same proce- 
dure as described in Q is followed. For every wounded 
nucleon, a contribution to the energy density with Gaus- 
sian shape (in x and y) and width oo = 0.3 fm is added. 
The amplitude of the Gaussian was adjusted to yield the 
same average multiplicity distribution as in the case with 
average initial conditions. We tested the effect of partial 
scaling with binary collisions and found that it has very 
little effect on the flow observables in the event-by-event 
calculation. The reason for this is the appearance of "hot 
spots" that work against the increase in initial eccen- 
tricity, because their expansion is not necessarily aligned 
with the event-plane. 

The prescription described in Refs. [23T - [29j is used to 
initialize the longitudinal profile, both for average initial 
conditions and for event-by-event simulations. It consists 
of two parts, namely a flat region around rj s = and half 
a Gaussian in the forward and backward direction: 



H(r) s ) = exp 



2** 



\Vs\ ~ ??flat/2) 



(12) 



The full energy density distribution is then given by 

s(x, y, n s ,b) = e H(n s ) n WN (x, y, 6)/n WN (0, 0, 0) . (13) 

For the LHC scenario, the parameters ijant and a v are 
respectively set to 10 and 0.5, in order to reproduce pre- 
dictions in [30j | . All parameters for RHIC are the same 
as in parameter set Au-Au-1 in 

The equation of state used in this work is that of the 
parametrization "s95p-vl" from Ref. [3l[, obtained from 
interpolating between lattice data and a hadron reso- 
nance gas. 

A Cooper-Frye freeze-out is performed, using 



E 



dN 



dN 



d 3 p dypxdpTC 



.9, / /Kp>^% , (14) 



where gi is the degeneracy of particle species i, and £ 
the freeze-out hyper-surface. The distribution function 
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FIG. 1. (Color online) Charged hadron multiplicity as a func- 
tion of the pseudo-rapidity n for most central (0-5%) and mid- 
central (10-20%) collisions, using Tfo = 136 MeV. 



is given by 
/(« M JV) = /o(«"Pm) 



1 



1 



(2tt)3 exp((u^ - /ii)/Tpo) ± 1 



(15) 

where pi is the chemical potential for particle species 
i and Tpo is the freeze-out temperature. In the finite 
viscosity case we include viscous corrections to the dis- 
tribution function, / = /q + Sf, with 



Sf = f (l ± /„)?>%, 



1 



2(e + V)T 2 



(16) 



where W is the viscous correction introduced in Eq. ([2]) . 
Note however that the choice Sf ^ p 2 is not unique [32| . 

The algorithm used to determine the freeze-out surface 
£ has been presented in 0]. In the case with average 
initial conditions, we include all resonances up to 2 GeV, 
performing resonance decays using routines from 33 - 3(| 
that we generalized to three dimensions. For the event- 
by-event simulations we only include resonances up to 
the 0-meson. Wc have verified that the effect of including 
higher resonances on vi is negligible in the case of average 
initial conditions, such that neglecting them in the event- 
by-event case is not expected to change the results for 
flow observables. The thermalization time used for the 
conditions at the LHC is To = 0.3fm/c (compared to 
t = 0.4fm/c for RHIC). 

For the sake of reference, the pseudo-rapidity depen- 
dence of the charged particle multiplicity distribution is 
shown in Fig.[T] The spectrum was obtained by fitting 
to predictions in [3p| and normalizing the multiplicity at 
r/ = to experimental results from [37 1. 

In the event-by-event simulation the same method as 
in is applied, where the flow coefficients 



v n = (cos(n(0 - ip n ))) 



(17) 
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FIG. 2. (Color online) Charged hadron V2 in 10-20% central 
collisions for average initial conditions (avg) and event-by- 
event simulations (e-b-e). The bands indicate the statistical 
error. Experimental data from ALICE [lj]. 
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FIG. 3. (Color online) Charged hadron V2 in 30-40% central 
collisions for average initial conditions (avg) and event-by- 
event simulations (e-b-e). The bands indicate the statistical 
error. Experimental data from ALICE [lj]. 



are measured with respect to the event plane, denned by 
the angle 



ib - 1 arctan ^ sin ( n ^> 
Yn — aician - - — 

n {pt cos(n<p)) 



(18) 



The weight pr is chosen for best accuracy [38( . 

In Figs.[2]and[3]wc present the elliptic flow V2 as a func- 
tion of transverse momentum px obtained in the event- 
by-event simulation, compared to the average initial con- 
dition case and data from ALICE [l|. 

The used freeze-out temperature Tpo = 135 MeV is 
the same as that used to produce Fig.Q] Event-by-event 
fluctuations have opposite effects in the two different cen- 
trality bins, increasing vi in more central collisions, but 
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FIG. 4. (Color online) Charged hadron V2 from event-by- 
event simulations as a function of centrality. 



decreasing it in 30-40% central collisions. This has been 
observed and explained previously at RHIC energies in 
Q . Overall, event-by-event fluctuations and finite viscos- 
ity improve the agreement with the experimental data 
for the 30-40% central bin. For p T < 0.7 GeV the ex- 
perimental data is underestimated. A similar mismatch 
to experimental data was found in [III 13 1, where [l3| 



explains the low pt V2 by including non-thcrmalizcd par- 
ticles from jet fragmentation. For 10-20% central colli- 
sion the underestimation of vi at low px is even larger. 
Note that we use the event-plane method to determine 
V2 , while the ALICE data is obtained using the four par- 
ticle cumulant method. This method eliminates most 
of the non-flow contributions but also minimizes effects 
of fluctuations [39]] . Using the event- plane method or 
the two-particle cumulant leads to about a 10% larger 
V2 [!(|. Potentially, initial conditions using the Monte- 
Carlo KLN Model [ill E3 > which increase the effective 
initial eccentricity, can improve the agreement with ex- 
perimental data by leading to larger V2- 

Fig. 2] shows V2 as a function of centrality compared to 
two particle cumulant f2{2} and four particle cumulant 
v 2 {4} from ALICE The results reflect the slighlty 
too low V2 at low pt seen in Figs. [2] and [3] 

The flow coefficient V2 for average initial conditions 
and a freeze-out temperature of Tpo = 160 MeV - just 
below the range of the cross-over transition - is shown 
in Fig. [5] It is very similar to the result obtained with 
Tpo = 137 MeV, which demonstrates that the elliptic 
flow is completely built up already at the earlier time, ex- 
plaining the small difference between the V2 (pt ) at RHIC 
and LHC. 

To make this point more clear, the time-dependent ec- 
centricity 
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FIG. 5. (Color online) Charged hadron V2 as function of 
transverse momentum pr for the final freeze-out temperature 
Tfo = 136 MeV and a freeze-out at the temperature below 
the cross-over region Tfo = 160 MeV. The difference is mini- 
mal, indicating that all the elliptic flow builds up during the 
QGP-phase. 



and momentum anisotropy of the system at midrapidity 



_ (T x 



Tyy) 



(20) 



are shown in Fig. [51 This illustrates how the momen- 
tum anisotropy is almost entirely built up after 5 fm/c at 
both RHIC and LHC and the (slightly) longer lifetime at 
LHC does not play a role for w 2 . We also find that while 
the system becomes just isotropic at freeze-out at RHIC, 
the eccentricity changes sign approximately 5 fm/c before 
freeze-out at LHC. Turning to the temperature evolution 
in the center of the system in Fig. [71 one concludes that 
at the LHC, almost all elliptic flow is built up in the QGP 
phase. The figure also shows that the lifetime of the QGP 
phase is approximately 40% longer at the LHC. The in- 
clusion of viscosity does not have a large effect on the 
temperature evolution at the center of the LHC fireball. 

These findings arc in line with earlier calculations using 
2+1D ideal hydrodynamics [43|, where a small decrease 
of «2 as a function of px was predicted, when going to 
higher initial energy densities. 

Finally, we present predictions for the triangular flow 
coefficient v^(pt) for LHC energies in Figs. [8] and [9l The 
results for rj/s = 0.08 are again remarkably similar to 
those at RHIC energies [H, while the difference between 
curves computed with different r/fs is smaller than at 
RHIC energies. 

We have presented first event-by-event calculations of 
elliptic and triangular flow at LHC energies of y/s = 
2.76 TeV. The data seems to be best described with 
a viscosity to entropy density ratio of rj/s = 0.08 or 
smaller. A value of rj/s = 0.08 was also found to describe 
RHIC data best earlier Generally, the Monte-Carlo 
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FIG. 6. (Color online) Time evolution of the eccentricity 
and momentum anisotropy in 10-20% central collisions at 
RHIC and the LHC. In both cases the maximal momen- 
tum anisotropy is built up almost entirely during the first 
5 — 6 fm/c. At LHC the spatial anisotropy drops even slightly 
faster. This demonstrates why the elliptic flow as a function 
of pr is very similar in both experiments. 
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FIG. 7. (Color online) Time evolution of the temperature in 
the center at RHIC and LHC. The horizontal line indicates a 
temperature of 195 GeV, which is in the center of the cross- 
over region defined in [311 ] . 



Glauber initial conditions tend to lead to relatively small 
i>2 in the 3+1D event-by-event simulation compared to 
experimental data from ALICE It was found that, 
in spite of the longer lifetime of the system at LHC, the 
anisotropic flow is not much larger than at RHIC: it is 
built up almost entirely within the first 5 — 6 fm/c, and 
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FIG. 8. (Color online) Charged hadron v-j, in 10-20% central 
collisions. 
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FIG. 9. (Color online) Charged hadron «3 in 30-40% central 
collisions. 



at the LHC the system even acquires a negative eccen- 
tricity for the last 5fm/c of the evolution, reducing the 
momentum anisotropy. The low value of r\js extracted 
from RHIC has been interpreted as the formation of a 
"perfect liquid"; it would seem this description is also 
appropriate at the LHC. A detailed combined compari- 
son of experimental spectra and higher harmonics using 
the techniques presented in this work will be able to de- 
termine the initial conditions and their fluctuations as 
well as the shear viscosity to entropy density ratio. 
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